Coupled S = I Heisenberg antiferromagnetic chains in an effective staggered field 
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We present a systematic study of coupled 5=1/2 Heisenberg antiferromagnetic chains in an 
effective staggered field. We investigate several effects of the staggered field in the higher [two or 
three) dimensional spin system analytically. In particular, in the case where the staggered field 
and the inter-chain interaction compete with each other, we predict, using mean-field theory, a 
characteristic phase transition. The spin-wave theory predicts that the behavior of the gaps induced 
by the staggered field is different between the competitive case and the non-competitive case. When 
the inter-chain interactions are sufficiently weak, we can improve the mean-field phase diagram by 
using chain mean-field theory and the analytical results of field theories. The ordered phase region 
predicted by the chain mean-field theory is substantially smaller than that by the mean-field theory. 

PACS numbers: 75.10.Jm 
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I. INTRODUCTION 

The effects induced by magnetic fields in magnets have 
been a subject of theoretical research interest for a long 
time. In particular, recently the magnetization processes 
of various spin chains and ladders have been investigated 
intensively. Owing to the progress of the various experi- 
mental methods, there has been an increasing connection 
between the theories and the experiments. In such a con- 
text, one of the new attractive subjects in magnetism is 
the effects of a staggered magnetic field, namely, a mag- 
netic field which changes direction alternatingly. While 
it may sound unrealistic, there exist at least three mech- 
anisms generating the staggered fields in real magnets, 
as discussed in Refs. U andQ- 

The first mechanism is due to the staggered gyromag- 
netic {g) tensor, which can be present if the crystal struc- 
ture is not translationally invariant. The staggered g ten- 
sor g^^p is defined in the coupling between the spin and 
the external magnetic field (Zeeman term) as 



(1) 



where H =^ (H^, Hy, Hz) is an applied uniform magnetic 
field and Sj is the spin operator of the local magnetic mo- 
ment. Here Hag^'^p is nothing but an effective staggered 
field. In addition, the staggered field may also arise from 
the staggered Dzyaloshinskii-Moriya (DM) interaction'^'^ 



X S. 



(2) 



which can be present if the crystal symmetry is suffi- 
ciently low. It is shown in Refs. Uind|3 that in the pres- 
ence of the staggered DM interaction along the chain, an 
applied uniform field H also generates an effective stag- 
gered field h oc D X H . Several quasi-one-dimensional 
Heisenberg antiferromagnets are now known to have the 
staggered field due to the above mechanisms. The well 



known examples are Cu-benzoate,^'^'* [PM- Cu(N03)2 • 
(H20)2]„ (PM=pyrimidine)^ and Yb4As33ii^ All of 
these have low-symmetry crystal structures which allow 
a staggered g tensor and a DM in teraction along the 
chain.ii It is expected in Refs. and Isl that the stag- 
gered field induces an excitation gap in the S — 1/2 
Heisenberg antiferromagnetic (HAF) chain, which should 
be otherwise gapless. The excitation gap caused by the 
staggered field is indeed found in these materialsiSiSiiS 
Moreover, low-temperature anomalies in physical quan- 
tities, such as the susceptibility and the electron spin 
resonance line-widthjiS4i& are also successfully explained 
as effects of the staggered field^. Thus it is confirmed 
that they are described by an 5 = 1/2 HAF model with 
an effective staggered field. 

There is another, rather different, mechanism to gener- 
ate a staggered field. Let us suppose that the system con- 
sists of two sublattices, with a weak inter-lattice coupling 
and strong intra-lattice one. If one of the sublattices is 
Neel ordered, the inter-lattice coupling, as a mean field, 
could give an effective staggered field on the other sub- 
lattice. The realization of this scenario is in R2BaNi03 
where R is a magnetic rare earth, and the R-ion lattice 
provides a staggered field for Ni chains {S = l)iiLi& 

Actually all the materials discussed above are highly 
one-dimensional (ID). However, at lower temperature 
and lower energy, the inter-chain interaction will eventu- 
ally be dominant. In addition, there are reports on a few 
materials [CuCl2 • 2DMS0 (DMSO=dimethylsulphoxide) 
(Refs. .19.20.21.1 and BaCu2(Sii_xGex)207 (Refs. Q^] 
which seem to have an effective staggered field and also 
a relatively large inter-chain interaction. Therefore the 
work including the inter-chain interaction could be rele- 
vant for experiments. 

Given these backgrounds, in the present paper, we 
would like to clarify the characteristic roles of staggered 
fields in higher dimensional spin systems. In this pa- 
per, we are concerned with dimensions higher than 1 but 
still realistic in condensed matter physics, namely, two or 
three dimensions. However, most of the analyses in this 
paper apply straightforwardly even to four dimensions or 
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higher. 

Varieties of spin models with effective staggered fields 
are conceivable. As a simplest model including both the 
staggered field and the inter-chain coupling, in this pa- 
per, we concentrate on the following 5=1/2 spatially 
anisotropic Heisenberg Hamiltonian 



H 



H = 



where Sij^k is the spin f/2 operator on r = (i, j, k) site. 
The coupling constants are restricted to J > \J±\ > \J'j_\, 
and thus the i direction is the strongly antiferromagnetic 
(AF) coupled one. The system with = is 2D, in 
which the index k vanishes. The last two terms represent 
the uniform and staggered Zeeman terms respectively. 

In our model Q, one can immediately find that when 
the inter-chain interactions are AF, they compete with 
the staggered Zeeman energy, while in the ferromagnetic 
(FM) case, both the interactions and the staggered field 
h jointly make a Neel state stable. Let us refer the AF 
case as the competitive case, and the FM case as the 
non-competitive one. As will be explained later on, we 
predict that the competition brings a second-order phase 
transition in the competitive case. Its emergence is one 
of the most characteristic effects of the staggered field in 
our higher dimensional spin model. 

The rest of this paper is organized as follows. In Sec. ITU 
we apply mean-field theory to the model Q . In the com- 
petitive case a phase transition is predicted. Since we 
are primarily interested in the transition, which is char- 
acteristic for the higher dimensional system, in the later 
sections we will mainly discuss the competitive case. The 
non-competitive case is touched as a comparison to the 
competitive case. Besides the phase diagrams, the mean- 
field magnetization curves and critical exponents are de- 
rived from the self-consistent equations. In Sec. lIIII using 
linear spin-wave approximation, we derive the spin-wave 
dispersions in the the competitive case and in the non- 
competitive case. As a result, we find that the excitation 
gap induced by the staggered field behaves differently be- 
tween the competitive case and the non-competitive case. 

In Sec. IIVI we improve the mean-field phase diagrams 
by using chain mean-field theory.^'^'^'* The latter is ex- 
pected to be superior, when the inter-chain interactions 
are weak and the effective ID model can be solved ex- 
actly. The improved diagram shows that in the weakly 
coupling region of the competitive case, the ordered phase 
becomes much narrower than the mean-field prediction. 
In the last section, we summarize those results and dis- 
cuss future problems. In Appendix, the details of the 
spin-wave results are given. 
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FIG. 1: Directions of the magnetic fields and the spin mo- 
ments in the 2D competitive case. The short black arrows are 
the spin moments projected onto spin xy plane in the ordered 
phase (spontaneous symmetry breaking phase, see text) ex- 
pected by the MFT. The gray arrows indicate the direction 
of the staggered field h. The uniform field H is applied per- 
pendicular to this paper. 



II. MEAN-FIELD THEORY APPROACH 

In this section, we treat the model Q within mean- 
field theory (MFT) framework. We first discuss the com- 
petitive case, and then touch the non-competitive case 
briefly. 

In the competitive case, considering the advantage of 
both the inter-chain energy and the staggered Zeeman 
energy as well as the intra-chain coupling J, we can ex- 
pect that the spin moment turns as Fig. ^ a-t sufficiently 
low temperature and small fields comparable to the inter- 
chain couplings. The assumed spin moment is 



i,j,k/MFT 



( i-iym.,, i-iy+^+'^my, m, ) . (4) 



The choice of the mean field Q is presumably valid as 
long as J is sufficiently larger than J± and J'^. On the 
other hand, for instance, in the limiting case: J — > 0, 
where the model is reduced to a 2D or ID AF one with 
a uniform field, there are possibilities that {Sj'^) are in- 
homogeneous along the j or k directions, and thus Eq. 
Q is invalid. In this paper, the MFT in the competitive 
case is performed only within the mean field Q). 

The minimal condition for the mean-field free energy 
gives the following self-consistent equations: 
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tanh(/3e), 



m 



tanh^(/3e) 



(5a) 
(5b) 
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where 



e 



h/2 



{J + + J'Jmy 
-{J + J^ + J'^)m,+H/2 ' 

(4 + 4 + ^1)'^' 



(6) 



and P and m are, respectively, the inverse temperature 
and the total magnetization per site. The numeri- 
cal solutions of Eq. (|5bp are given in Fig. 13 They in- 
dicate that there is a second-order phase transition, and 
the corresponding order parameter is the y component of 
the spin moment niy = \{Sfj^.)\- Going back to Fig. ^ 
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FIG. 2: Magnetization curves of the competitive case for 
(J, J± + J'j^, H) ^ {1, 0.2, 0.2). These are obtained by the nu- 
merical iterative method for Eq. llSbl . The upper two parts 
(a) and (b) are in h = 0.05 and h = 0.35, respectively. The 
lower two parts (c) and (d) are in ksT = 0.3 and ksT = 0.7, 
respectively. 

one sees that the phase with finite ruy breaks the transla- 
tional symmetry in the weakly coupled direction. In the 
following, we call this phase as the SSB (spontaneous 
symmetry breaking) phase. The other phase, in which 
the spins are aligned to the field with my = 0, will be 
called as the symmetric phase. We emphasize that the 
transition between these two phases occurs only in the 
high dimensional spin systems and in the presence of the 
competition. From these results, we can illustrate the 
variation of the spin moment when the staggered field h 
is increased gradually at small T and H with Fig. |3| 
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FIG. 3: Variation of the spin configuration in a plaquette 
when the staggered field is increased gradually in the 2D com- 
petitive case. The black arrows represent spin directions pro- 
jected onto spin xy plane. The gray arrows are the staggered 
field h. 
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FIG. 4: Schematic mean-field phase diagram in the competi- 
tive case. 



In the SSB phase, the relations nix 

. _ H 



and 



hold within the MFT. Inserting these 
0, we obtain the 



'"2 - 4(.7+ Ji+J^) 

into Eq. (|5b|l and taking the limit rUy 
mean-field critical surface in the space {ksT^ iJ, h): 



H, 



2 = - tanh^ 



l3,{J + J^ + J'^)Jhl + H^\, (7) 



where hr 



Ha 



and the sub- 



script c represents critical values. It can be simplified in 
the cases T — 0, h = and H — 0, respectively, as 



i tanh 



he = h tanh |/3c(>/ 



f /9c He 

V 4 



Jl 



(8) 



Thus the mean-field phase diagram can be represented 
as Fig. 21 Using the critical condition, one can calculate 
some critical exponents within the MFT. Near the critical 
surface in the SSB side, the order parameter, the off- 
diagonal uniform and staggered susceptibilities: Xu = 



4 



and Xs = behave , respectively, as rUy ~ {Ac — 
A)'^, Xu --{Ac- A)-! and Xs --{Ac- A)-<' where A 
stands for T, H or /i. The critical exponent (3 is found to 
be the conventional mean-field value 1/2. On the other 
hand, both 7 and 7' turns out to be 1/2, which is different 
from the standard MFT result 1. This is because ray is 
perpendicular to H and ft,, and thus the latter are not the 
conjugate field as in the standard case. The mean-field 
energy per site is given as 



CMFT = ^{^( 

-f(J_L + J'A_){ml - ml 



l + niy — ml) + Hmz + hmx 

(9) 



From this, one can easily confirm that the critical expo- 
nent for the specific heat is zero. 

Now, wc turn to the non-competitive case. Because 
of no competitions, no singular phenomena occur when 
h ^ 0. Canting of the spins in zx plane lowers both the 
inter-chain interactions and the Zeeman energies. Thus 
the expectation value of the spin moments can be put as 



(<^ij,fc)MFT = ( {-^Ym^, 0, m^ ) 



(10) 



The MFT in this case gives the self-consistent equations 
^^(.-) = «(^)/2e')tanh(/3e') where e^^^^ = [{-)J+\J^\ + 
\J'J]mx(,) + h{H)/2 and e' = {el + tlfl'^. At /i = 0, the 
system reduces to a conventional AF magnet in a uniform 
magnetic field. Hence, there must be a phase transition 
which divides the AF and paramagnetic phases, charac- 
terized by the order parameter m^- The critical line is 
given by 



|^ = itanhU(J-H|J^| + |Ji|)|^ 



(11) 



The phase diagram and the variation of the spin moment 
are drawn in Fig. In the three-dimensional parameter 
space {ksT, H, h), the AF phase gives a first-order phase 
transition plane. 



III. LINEAR SPIN-WAVE APPROXIMATION 

IN T = 

With the MFT described in the preceding section, we 
investigate the effects of the quantum fluctuations in both 
the competitive and the non-competitive cases, at T = 0. 
The standard linear spin-wave approximation, based on 
the Holstein-Primakoff transformation (HPT) (Ref. f25l) 
is employed. The detailed results are given in Appendix 

El 

First we discuss the SSB phase of the competitive case, 
which is the main subject. In the HPT, we replace the 
spin operator with a boson annihilation (creation) oper- 
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FIG. 5: (a) Schematic mean-field phase diagram of the non- 
competitive case in h — 0. (b) the variation of the spin 
moment when the uniform field is increased gradually in 
ft = T = 0. 



ator c (c^) as follows: 



S'hp — 



S 



|(ct 



(12) 

where S is the spin quantum number generalized from 
S — 1/2 and the sign « denotes the leading approxima- 
tion in the 1/5* expansion. The above HPT is useful if 
the spin points to the y direction in the classical ground 
state, as the bosons then represent quantum fluctuations. 

In the SSB phase, actually, a canting structure Q is 
expected in the classical ground state, which is equivalent 
to the MFT at T = 0. The canted spin moments may be 
expressed by the angles {9, (f>) as 



(—1)' cos 9 sin (f> 
(S^jj,fc)=m( (-1)*+^+'^' cos 61 cost 
sin 6* 



(13) 



Thus in order to apply the HPT to the present case, we 
use the representation 



r.3,k R.{{-iy+' + ''9) X 

Rz[{-ly+''{-4>)+S^+J+kMd ^)5hp,(14) 

where the operator Ra{P) represents a rotation about a 
axis by angle /3. The operation (|14|) is described in Fig.El 
For the original spin model (Q, these rotations corre- 
spond to a unitary transformation. After these trans- 
formations, the leading order terms in l/S* is retained to 
give a solvable Hamiltonian which is quadratic in bosonic 
operators. 

The canting angles of the classical ground state are 
given by minimizing the classical Hamiltonian, which is 
equivalent to the MFT energy, Eq. (|5Jl. The canting an- 
gles in the classical ground state are thus given as 



COS( 

sin 6 



ci sm 9ci = 



H 



4S(J+,/_L + ,/^ 



=1 - 4S(Ji + ,/' )• 



(15) 
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FIG. 6: Definitions of the canting parameter. In the SSB 
phase, as in (a), each spin direction is given by two rota- 
tions (1) -Rz(— 0) and (2) Rx{0) applied in this order to spin 
pointing to y direction. In the symmetric phase or in the non- 
competitive case, spins are canted in the zx plane as shown 
in (b). 




FIG. 7: Magnon dispersions of the SSB phase in the 2D com- 



The gray circles represent the gapless points. 



For the resulting quadratic Hamiltonian iJnp, we petitive case for (S, J, Jx) = (1/2, 1, 0.1) and 
perform a Fourier transformation (FT) '^i j k ~ 
A^^*^/^ e''^ ''ct, where N is the linear system size, d 

is the dimension of the system, fc^ (|fca| < T^/aa) and 
are, respectively, the wave-number and the lattice con- 
stant for the a-direction. A four-mode Bogoliubov trans- 
formation (BT))2Si which mixes cj, c_g, cj_-, c and 
their Hermitian conjugates (tt = (0, Tr/a^, Tr/a^)), leads 
to the diagonalized form 



k 



uj{k) clc^ + const. 



(16) 



with a single band in the first Brillouin zone, where Cj: 
is the magnon annihilation operator and the direction 
corresponds to the strongly AF one. The explicit results 
on the dispersion uj{k) are lengthy and thus are given in 
Appendix ^ Here we discuss physical implications of 
our results. 

First, the obtained dispersion satisfies uj{k) > for all 
values of parameters. This implies that the SSB phase, 
which appears as the classical ground state, is stable 
against quantum fluctuations, at least in the lowest order 
of 1/5. Some representatives of the dispersion arc shown 
in Fig. [3 At zero field {H ^ h ~ 0), there are two gap- 
less points kfi = (0,0,0) and kn = (Tr/a^;, Tr/oy, tt/o^), 
with linear dispersions in the neighborhoods. Let us de- 
fine A;i = Lo{kfi) and Ah = Lo{kH)- Since our results in 
Appendix lA II indicates that Ah (Ah) is non-vanishing 
only when h {H ^ 0), we call it as /i(iJ)-induced 
gap. The true excitation gap, namely, the minimum ex- 
citation energy, is given by A = min(A/i, A//). Thus 
the gap A vanishes exactly as long as either h or H 
remains zero. This is contrasting to the S — 1/2 ID 
HAF model, where the staggered field alone induces the 
gap.^^ The gapless excitations are identified as Nambu- 
Goldstone (NG) modes. Indeed, when either H oi h is 
zero, the Hamiltonian has a continuous U(l) symmetry, 
which is broken spontaneously in the SSB phase. 

In the non-competitive case, the magnon dispersion 
relation (D(fc) is given in Eq. l|J20|) in AppendixlX2l Now 



the staggered field alone can open the gap, because the 
ground state does not break any continuous symmetry 
spontaneously. On the other hand, the system remains 
gapless at fc = due to the NG mechanism ii h — and 
H is not too large. The /i-induced gap is thus defined as 
A^ =LD(fc = 0). 

While the /i-induced gap is defined for both cases, there 
is a characteristic difference in the h dependence of the 
gaps. In the limit of small h, Ah ~ h^^'^ for the non- 
competitive case, but Ah h for the competitive case. 
In other words, one can say that the /i-induced gap in 
the competitive case opens more slowly than in the non- 
competitive case. This is naturally understood because 
the competition between J± (or J^) and h weakens the 
effect of the external symmetry breaking hy h. It also 
implies that the ground state is more stable in the non- 
competitive case, against quantum and thermal fiuctu- 
ations. The opening gaps are drawn in Fig. |S1 On the 
other hand, in the limit of small H, A^ ^ H. Similarly 
to the case of Ah, this may be interpreted as a result of 
the competition between the uniform field and the AF 
couplings. 

We expect the spin-wave theory for the gaps to be 
qualitatively correct even for S = 1/2. As discussed in 
Appendix I A 21 the spin-wave dispersion of the symmetric 
phase can be obtained by the replacement (| J^|, | J^|) ^ 
(— J^,— ,/^) in the dispersion of the non-competitive 
case. 

Finally let us discuss the relation of the present results 
to the spin-wave theory in the ID model {J\ = J'j^ — 0) 
with the staggered field h discussed in Ref. Q. The spin- 
wave dispersion a)(fc) for the non-competitive case does 
approach smoothly to the ID dispersion [Eq. (3.11) in 
Ref. 5], when J_l, 0. On the other hand, the dis- 

persion uj{k) for the competitive case apparently does not 
reduce to the ID one in the limit of J_l, 0. This 

is because uj^k) is the dispersion of the magnon excita- 
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FIG. 8: Ah and in the 2D case for {S,J,J±,H) = 
(1/2,1,0.1,0). In the present case where H = 0, the gaps 
have simple forms: Aj, = (1 + J/ J±)^^^h and Aj, = [1 + 
h/{ASJ)]^/'^VISJh (see Appendix). 



tion in the SSB phase, which is absent in the ID model. 
In fact, for any finite h, the SSB phase is realized [and 
hence uj{k) is applicable] only when J± and are above 
the critical values. Thus by decreasing Ji_{J'^_) at a fixed 
h, the system undergoes a phase transition into the sym- 
metric phase. At the transition the dispersion should also 
change drastically. In the symmetric phase, similarly to 
the non-competitive case, the dispersion does approach 
continuously to the ID dispersion. 



3,k 



-{-iy+^[h' + 2(j^ + J'^)my]{-iysl^^ 

-[h-2{J^ + J'^)m^]{-iySl^^^ 
-[i7-2(J^ + Ji)m,]5f^^.„ (17) 

where we introduced an infinitesimal staggered field 
{-\y+i+^h' parallel to the order parameter. The CMFT 
requires that the mean fields mx,y,z are equivalent to the 
corresponding moments of the effective chain Hj^k- It has 
the two effective staggered fields = h — 2{J± + J'jjmx 
and hy = h' + 2{J± + J'^)my as well as the effec- 
tive uniform field ^ H — 2{J± + J'^JrUz- Clearly 
it is sufficient to consider one chain where j + k is 
even, and we represent it as Hid. Within the linear- 
response theory, the moment | (5'^^) | , in which (• ■ ■ ) stands 
for the mean value of Hiu, can be approximated by 
xl^{Hz,hx, 0)hy where the staggered susceptibility is de- 



fined as Xy {Hz,K,hy) E 

requirement leads to my 
to 



d\{Sl)\ 



Therefore the above 



Xy^hy which is transformed 



IV. CHAIN MEAN-FIELD THEORY 
APPROACH 



l-2(Ji+,/' )xi 



h'. 



(18) 



In this section, we reconstruct the phase diagram of 
our model using chain mean-field theory (CMFT). In 
the CMFT, weak couplings among the chains are treated 
with a MFT, and the resulting effective ID problem is 
analyzed as precisely as possible. If the ID problem can 
be treated exactly, the CMFT is expected to be much 
reliable when the one dimensionality is strong enough as 
in the case of Cu-benzoate, since it includes the fluctu- 
ations in the strongly coupled direction correctly. The 
usefulness of the CMFT has been demonstrated in sev- 
eral applicationsi^^i^^i^^ 

In Sec. IIV Al we discuss how the CMFT determines 
the phase transition for our model lj2Jl. Sec. IIV 51 is a 
brief overview of susceptibilities of the S* = 1/2 HAF 
chain which are necessary to the CMFT. In Sec. IIV CI we 
present the CMFT phase diagrams and compare them to 
the MFT ones. 



A. CMFT for Our Model 

Let us derive the effective ID model for our system (O, 
within the CMFT. 

In the competitive case, we consider the symmetric 
phase side for convenience. The mean-field procedure for 
the weak inter-chain couplings replaces them with the 
effective external fields. Thus the resulting Hamiltonian 



Similarly and ruz can be determined by the CMFT 
as well. The critical condition of our model is 



l-2(Ji + Ji) xl''{Hz,hx,0) = 



(19) 



At this point, one sees that the original 3D or 2D suscep- 
tibilities and moments can be described by those of the 
effective chain. 

The staggered field in Hir, has both x and y com- 
ponents: hx and hy. By the infinitesimal rotation Rz{6) 
about z axis by angle 5, where sin (5 = —[hy/{hy + h'^y^^'^, 
the effective Hamiltonian is simplified as 



ID 



HzS^ 



Ki-^ys'r, (20) 



with the one-component staggered field h'^ — hx/ cosS. 
In order to obtain an alternative formula determining my 
instead of Eq. (|18|l , we focus on the relation 

(5f )' = cos (5 {S^)-sm6 (Sf) 

= (— l)*(cos (5 m^; — sin (5 JTij,), (21) 

where (•••)' represents the expectation value of H[y) 
and the second equality is caused by the self-consistency 
of the CMFT. Here we define the susceptibilities of 



the HAF chain ^ as x}P(T,Hz,h'J 



and 



7 



xl^'iT^Hz^K) ^ ^ where m^., = 



Within compactification radius of the effective boson field the- 



the linear-response theory Eq. 12111 is reduced to 

{T, Hz,h'^) h'^^ cos (5 - sin 5 my . (22) 

Let us recall that S is defined by and my, and that 
mx,z are determined by the CMFT. Consequently, my 
can be determined as a solution to Eq. 1)22(1 . The condi- 
tion that my diverges would determine the phase transi- 
tion. 

B. Susceptibilities of 5 = 1/2 AF chains 

In order to solve Eqs. (|19|l or H22|l in terms of 
my, we need the explicit forms of the susceptibili- 
ties: xlP and xl^ ■ Here we briefly summarize the 
known resultai^*22i22i^ on these quantities, obtained by 
the bosonization technique. 

In the absence of the staggered field h'^, the low-energy 
effective theory of of the Heisenberg chain (|20|l is given by 
a free boson field theory, and the uniform susceptibility 
at zero temperature is obtained as 



(23) 



(27r)2i?(i7)2«(iJ) ' 
where = H, and R and v, respectively, are the 



ory 



29,30,31 g^j^j ^]-^g spin-wave velocity. The exact values 
v{H) and R{H) as functions of the uniform field H are 
given by a solution of a set of the Bethe ansatz integral 
equations.'^— In the case oi H = Q, they are explicitly 
given as 



■Kja/2, i?=l/\/2^. 



(24) 



For a small uniform field = i?(<C J), the asymptotic 
behavior of the radius R follows^ 



27ri?^ 



1 - 



1 



2\n{J/Hy 



(25) 



The logarithmic temperature correction due to the 
mar gina l operator of the HAP chain is discussed in 
Ref. |33- Although the transverse staggered field h'^ in- 
duces a gap, it is expected to have little effectJ^ on the 
uniform susceptibility xV^ if ^'x small enough. 

Next we turn to the staggered susceptibility Xs^- In 
the case J h'^, ksT, the low-energy effective theory of 
the Hamiltonian H20|) is given by a quantum sine-Gordon 
field theory. Using the exact solutions of the HAF chain, 
the scaling arguments and Lukyanov-Zamolodchikov pre- 
diction24i2^, the staggered susceptibility Xs° for small 



Hz = H and h'=hatT = Q^-'^ 



IS given as 



L>[2(1 - 27ri?2)] 



-1/3 ttR^ 



-2(l-27ri?,^)/3 



iH = 0) 



(26) 



where the radius R is that of the model without the 
staggered field, and D = 0.3868 .... The first formula 
is actually valid for < 2 J (below saturation field), 
but the second is only so for _ff <C J. These formulas 
are correct in ksT h <^ J. (More precisely, within 
ksT <C Ah <C J where is the /i-induced gap in the 
HAF chain.^) 

In the intermediate temperature regime h, H <C 
ksT ^ J, where the temperature is larger than the in- 
duced gap, the staggered susceptibility may be approxi- 
mated by that for zero staggered field^ as 



x,i°(fci3r»/i) 



[ln( J/fc^T)] 
ksT 



(27) 



C. Phase diagrams in CMFT 



Employing the results of Secs. lIVXl and UVBl we study 
the phase diagrams, in particular for the competitive 
case, within the CMFT. Unfortunately, it can be applied 
only to several limited regions in the parameter space, 
where the susceptibilities of the ID model are obtained 
exactly. 

First, let us consider the region near {ksTc, 0,0). In 
the zero-field case, the effective model Hid has only an 
infinitesimal staggered mean field hy as long as it is in 
the symmetric phase. Therefore the formula (|27|) can 
give the self-consistent value of my, and Eq. H19|) is the 
critical condition which serves the critical temperature: 



fcsT, = 21? ( + Ji 
« 21? ( Jj_ + 



In 



In 



1/2 



.7 



2V (J^ + J'J 



1/2 



(28) 



where V = 0.2779 . 



which is reasonable when fcsTc, (Jj^ + J^) ^ J and 
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FIG. 9: (fcsTc, J± + J'^) in the competitive case for J = 1. 
(a) and (b) are, respectively, predicted within the CMFT and 
the MFT. It is confirmed that the SSB areas of the CMFT 
are considerably smaller than the MFT. 
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FIG. 10: {hc,J± + Ji) in the competitive case for J = 1. 
(a) and (b) are, respectively, predicted within the CMFT and 
the MFT. Similarly to Fig.^l the narrowing of the SSB areas 
occurs. 



agrees with the result of Ref. |3- Equation 128|) is also 
stable against a smaU staggered field h, because xl^ is 
independent of h in Eq. H27|l . These results are also 
valid for the non-competitive case, with the replacement 
J± + J'j_ I + I Jil- Figure El represents the compar- 
ison between the CMFT result (|28|l and the mean-field 
prediction: ksTc = {J + J± + J'^) /2. 

Next we investigate the neighborhood of (0,0, /ic). In 
this case {H = 0), Hid has two staggered fields. There- 
fore through the rotation of Sec. IIV Bl the order param- 
eter iriy can be fixed by Eq. 122() . Inserting the third 
formula of Eq. (|26|l into Eq. (|22|) . and performing the 
Taylor expansion of Eq. H2()|l around hy — 0, we obtain 
the linear-response relation: 



{|-l[ln(J/fe.)]-i}m, 



y /i-i{lO(,/i + J^)-[ln(J//i,)]~i}m 

Therefore the critical condition can be written as 



h'. (29) 



10(J±-K Ji)' 



(30) 



where we assumed he, {Jj_ + J'j_) ^ J. On the other 
hand, rrix should be fixed from the self-consistency of the 
CMFT as well. At T = 0, it leads to 



xl°{0,0,h-2{J± + J'Jm^) 



(31) 



Combining Eqs. H30|) and (|31|l . we obtain the critical stag- 
gered field: 



hr. = D' 



D' 



1/2 



2( Ji + Ji 



1/2 



c2(Jj_ + Ji 



1/2 



(32) 



where we kept only the leading order of {J± + J'^)/ J and 
£>' = (2 X lO^D^S-^)!/^ ~ 7.27- • • X lO^^. This is valid 
if he, (J± + Ji) ^ J, and is compared to the MFT result 
he = 2{Jj_ + Ji) extracted from Eq. © in Fig.UHl The 
CMFT correction to the MFT is found as a significant 
multiplicative factor D'[- ■ -J^/^. 

We consider furthermore the critical line in the T = Q 
plane in the limit of he <C He ^ J. In the sym- 
metric phase near this line, the effective model i/iD 
has all three kinds of the external fields. Hence rrix 
and TUz must be determined concurrently by the CMFT 
scheme. However in the present case he ^ He, may 
be estimated independently by taking an approximation 
xi°(0, H, h) ~ xi°(0, H, 0) which was given in Eq. O^ . 
From this approximation, rriz is also fixed as 



H 



^2 J(l - [ln(J/iJ)]-i) + 2( Ji + Ji) ' 



(33) 



which is justified in H, (Jj^ + Ji) ^ J. In the sufficiently 
small field case, J 3> i/(^ h), the logarithmic part can 
be dropped as well. From the second of Eq. H2t)|l and 
Eqs. H30|) . H31|) and H33|) . we obtain the critical line in 

T = 



he = D" 



J±+J\ 



, / 7r^/(l-[ln(.//g.)]"i)+2(./^+.7i) 



+ [\n{J/He)]-^ + 



1/2 



2{J± + Ji) 



(34) 



where D" = (lO'^D^e-^S-^Y^^ « 5.40 •• • x IQ-^. This 
condition is presumably suitable for fc^Ti ^ he <C He <C 
J. Thus, in order to obtain a more precise condition, it 
is necessary to adopt the first formula (|26|l and the exact 
values of R and v. Figure 1111 exhibits the comparison 
between this line and the mean-field prediction. 

Finally, we consider the point {0,He,0), where J/id 
has only the uniform field H^ — He — 2{J± + J'^)mz 
except for the infinitesimal field hy. It has been known 
from Bcthe ansatz that the magnetization saturates at 
the point iJ^ = 2J at T = in the HAF chain having 
only the field Hz. The transition point between the SSB 
and symmetric phase in this case should be identified 
with the saturation of the uniform magnetization. Hence, 
within the CMFT, the critical uniform field is given as 



He = 2J+(Jx+Ji). 



(35) 



The substitution J_i + Ji -(k-i| + ^il) gives the 
critical field of the AF-paramagnetic transition in the 
non-competitive case. 
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FIG. 11: Critical surface {he, He) in the competitive case for 
(J,T) = (1,0) obtained with the CMFT ^ (upper panel) 
and the MFT Q (lower panel). In the whole region of the 
CMFT panel, the critical region surrounded by the dashed 
ring {he ^ He <C J) is greatly reliable. 



From these results, we can compare the CMFT and the 
MFT phase diagrams. The comparison for the compet- 
itive case is summarized in Fig. 1121 In the case of weak 
inter-chain couphngs ( J ^ J±, J'±), the SSB phase of the 
CMFT is much smaller than one of the MFT. Especially 
there is a significant narrowing in ksT and h directions. 
As seen from Eq. (|18|) . this is because the phase transi- 
tion in the CMFT framework is driven by the divergence 
of the susceptibihty (in the present case, xl^) in the ef- 
fective chain, while temperature and the field (or rrix) 
suppress the divergence. On the other hand, the reduc- 
tion of the critical uniform field He is small in the weakly 
coupled case. This is because the uniform field competes 
with the intra-chain AF interaction as well as with the 
inter-chain interactions. 

Finally wc review the validity of the two theories. Since 
the strong one dimensionality is the basis of the CMFT 
procedure, it is expected that the CMFT is more reliable 
in the limit of J 3> J±tJj_- On the other hand, when 
J± and J'j^ are comparable to J, the special treatment 
of only one direction is unjustified. Therefore the MFT, 
which treats all couplings equally, is more reasonable for 
J ~ J±, J'j^. 



V. SUMMARY AND DISCUSSION 

We considered the effects of the staggered field h in 
an S* = 1/2 Heisenberg antiferromagnet Q in two or 
three dimensions. The system behaves quite differently 



ceo 

phase 



keT 



Critical surface in CMFT 



^ critical surface in MFT 

" symmetric phase 

FIG. 12: Schematic phase diagrams of the CMFT and the 
MFT in the competitive case. 



depending on whether the staggered field and the inter- 
chain couplings are competitive or not. In the com- 
petitive case, the appearance of a characteristic ordered 
(SSB) phase is predicted by the MFT. The SSB phase 
breaks the translational symmetry of the weakly coupled 
direction, and therefore it is peculiar to high dimensional 
systems. We also applied the CMFT to the model (O, 
and predicted that the region of the SSB phase becomes 
narrow in the CMFT scheme. The MFT and the CMFT 
are valid, respectively, in J J^, and in J ^ J ±,,7'^. 
The crossover behavior between these two regions can 
not be described by the mean-field type approachi^SiSSiSS 
It would require a more precise treatment of fluctuations. 

Moreover we studied spin- wave theory in both the com- 
petitive and non-competitive cases at T = 0. When 
the uniform field H is non-vanishing, the /i-induced gap 
opens as A/j ^ h in the SSB phase, while A;i ~ /i^/^ 
in the AF phase of the non-competitive case. This dif- 
ference reflects the partial cancellation of the staggered 
fleld effect due to the competitive inter-chain interaction 
in the SSB phase. The spin-wave dispersion in the SSB 
phase remains gapless due to NG mechanism even under 
a non- vanishing h. This is in contrast to the case of the 
ID model. 

Finally we comment on a few recent reports related 
to our study. In BaCu2Si207 reported in Ref. js^ both 
the staggered fleld and the inter-chain interactions are 
expected, as in our models. However, the effect of the 
exchange anisotropics, which is ignored in the present 
paper, is argued to be responsible for the observed two 
spin-flop transitions. Extending the present work to such 
a system would be an interesting problem in the future. 
Furthermore, in BaCu2(Sii_xGex)207rS the sign of the 
inter-chain interaction seems to depend on the doping 
parameter x. Thus, it could provide a realization of the 
competitive and non-competitive cases. 

Wang et al^ investigated an S = 1/2 AF ladder sys- 
tem with a staggered field. They argue that the competi- 
tion between the staggered field and the rung interaction 
brings a quantum criticality. It might be interesting to 



10 



compare our analysis on the higher dimensional system 
with theirs. 
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where * and T, respectively, stand for the complex con- 
jugate of each matrix component and the transpose of 
matrices, and the 4x4 matrices rjj: and the column 
four- vector are given as 



APPENDIX A: DETAILS OF SPIN- WAVE 
RESULTS IN SEC.irm 

Here we supplement the spin-wave results omitted in 
Seclml 



1. The competitive case 

We write down the details of the competitive case. Af- 
ter the FT of the boson operator Qj^fc, the spin- wave 
Hamiltonian can be expressed as the following matrix 
form: 



4^ = 



Ei(k) iEi{ky^,) 

-lEiiky,,) Fi{k) 



V -iEi{ky^,) 

f E2{k) iE3ik,)\ 

E2{k) iE^{k^) 

iE^{k,) F2{k) 

KiEsik,) F2(fc) / 

( ^—k ^k—T? ^—k — TT ) 



Fi{k) I 



(A2) 



4 



-AEi(k) 
(Al) 



Here -Ei,2,3.4 are defined as 



Ei{k) 

E2{k) 

Eaik^) 
Ei{ky^z) 



9cl COS^ 0cl cos kxCLj 



25*7(1 — 2 sin^ 9ci cos^ (f>ci) + 2S{J± + ./1)(2 cos^ flci cos^ (f>^i - 1) - 2SJ sin^ 9ci cos 

|-2S'(cos^ 6c\{l + sin^ 0ci) — 1)(>^± cos kyGy + J'^ cosfc^a^) -t- i7 sin ^ci cos 0ci + hshi(j)c\ 

SJ{cos^ 9ci — sin^ 6*01 sin^ (pci) cos kxttx + S cos^ 6ci cos^ 4'ci{J± cos kyUy + J'^ cos fc^a^) 
— S'Jsin 20ci sin0ci cos /c^^Oa; /4, 
S sin 26'ci sin (/)ci( Jj_ cos kytty + J'^ cos kzaz)/'^, 



/4, 



(A3) 



r 



and i^i,2(fc) = Ei^2{k — tt). Thus is Hermitian and -q^ 
is symmetric. Let us suppose that a four-mode BT 



(A4) 



where C-j^ 



[cr,c r,cr r -)"^ is a set of new bo- 

son (magnon) operators and AfBT(^) is an 8 x 8 matrix, 
diagonalizes the Hamiltonian as follows: 



ffHP 



fe 



^k 

n^. 



4 



4^1 (A:) 
(A5) 



where £7^ = diag(wi(A:),W2(fc),<j-'3(fc),CL'4(fc)). According 

to Ref. [2^ determining il,^ and Mbt(A:) is equivalent to 
solving an eigenvalue problem 



^k 



Ml^ = ( ) • (A6) 



The eight eigenvalues ztujj{k) [j = 1, 2, 3, 4) are given by 
±A/ {I = 1,2,3,4), where 



Ai(^) = [l/2G,--l/2(G|-4Gt)i/2]i/2, 

A2(^) = [l/2G,--l/2(G|-4G-)i/2]i/2, 

A3(^) = [l/2G,-+l/2(G|-4G+)i/2]i/2, 

A4(^) = [l/2G,~+l/2(G|~4G-)i/2]i/2^ (^7) 
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and Gr and are defined as 

^ k 



El 



El 



2El 



2El 



Gl = {El -El) 



{El - El) {Fl ~ Fi) 



-2{EiFi-E2F2){El+El) 
=p4 {E1F2 ~ E2F1) E3E4. 



(A8) 



The physical dispersion a;(/c) is given by either 4[Ai(A:) + 
X2{k)] or 4[A3(A;) + A4(fc))], depending on the value of k. 
In the vicinity of the point k — (fc^f), which is the 



gapless points when h = {H = 0), we find 



c^(fc) = 4[Ai(fc) + A2(fc)]. 



(A9) 



If only one of the external fields {H or h) is non- vanishing, 
the result is considerably simplified because we have 
E^ = E4 = 0. The spin-wave Hamiltonian can be ac- 
tually diagonalized by a simpler two-mode BT which is 
the same type as we need in the non-competitive case in 
Appendix lA 21 The resulting magnon dispersion is 



J 



ioik) = (2SJ{1 - 2sin^ 6*01 cos^ 0ci) + 2S{J±_ + J;^)(2cos^ ^cicos^ ^ci - 1) - 25Jsin^ ^cicos^ (t)ciCOS k^a^ 

-|-2S'(cos^ 6'ci(l + sin^ ^d) — 1){J± cos kyUy + J'^ coskzUz) + E j 

/ \ 2-1 1/2 

— ( 25J(cos^ 9ci — sin^ 9^ sin^ (f>ci) cos k^a^ + 2S'cos^ 9ci cos^ (t>ci{J± cos kyUy + J'^ cos kzaz)) ,(A10) 



r 



where T, = H sin 9ci cos 0ci when h — 0, and Y. ~ h sin (j)ci 
when H = Q. Of course, in these special cases, Eq. (|A9p 
reduces to Eq. (|A10|) near k — kh or kn- 

From the dispersion (|A9|) . let us estimate how (Ah) 
grows when a small h {H) is applied. To estimate Ah, it 
is sufficient to know the coefficients of Taylor expansion 
of i?i,2,3,4(^ = 0), G^^^ and - around h — 0. As a 
result, in the limit of small h, the gap behaves as 

A,, « 2^(J+Ji+Ji)l/2(./^ +.7^)1/2(1 _ij2)l/2 

-1/2 



1 _ 4^2 _ iJ-j^-j\f 

.j(.j±+j'j ^ 



where H — 



H 



and h — 



, (All) 



At H = 



4:S{,J+J±+J'j_) """"" ~ 4:SiJ±+J'j_}- 

0, Eq. HA11|) is reduced to the exact result A^ [1 + 
J/ {J± + J'j_)]^/'^h which is derived from Eq. IjAlOp and is 
drawn in Fig.|Hl Similarly to A^i, Ah can be estimated. 
The result is 



25(J+(Ji + Ji)(l-/i2))i/2 

xij^ + ,r^)'/\i-hY'^' 

{J-J_L-J'_^f 



1/2 



xH 



1/2 



(A12) 



which is reduced to the conventional result Ah — H 
when h = 0. Both gaps have linear field dependence. The 



results (|A11|) and (jA12p indicate that Ah {Ah) is non- 
vanishing only when /i 7^ {H 7^ 0). As a consequence, 
the true gap A = min(A;i, Ah) is zero when either of h 
or H is zero. 



2. The non-competitive case 

Here we summarize the spin-wave approximation on 
the non-competitive case. We find that it can be straight- 
forwardly applied also to the symmetric phase in the com- 
petitive case. 

According to the MFT and Fig. |H1 the spin configura- 
tion in the classical ground state can be written as 

52„j,fe - SRy{-9)x, (A13) 
S2n+i,3.k ~ SR,{7:)Ryi-9)x, (A14) 

where x is the unit vector pointing to x direction and 9 
is the canting angle parameter. The minimization of the 
classical energy determines 9 = 9c\ as 



cos 9c\ sm 9c\ = 



H 

4SJ 



cos ticl 



45J 



sm Oci- 



(A15) 



A standard spin-wave theory on this classical ground 
state gives the quadratic Hamiltonian in terms of bosons 



HP2 



E 



2 k 




2 k 



-A^ -I- const. 



(A16) 
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where 



- 2SJ cos29ci ~ 2SJ sin 6ci cos k^a^ 
+25*1 J_l|(1 — cos kyGy) 
+2S'|Ji|(l -cosfc^aj (A17) 
+H sin ^ci + h cos ^ci, 
Bk, = B-k, = ~S J cos^ 6ci cos kxa-x- 



Now we apply the two-mode BT 




(A18) 



where ui — vi = 1, to Eq. (jA16|) . The HamUtonian is 
then diagonahzed as 

HiiP2 — ^ (^{k) at dj: + const. , (A19) 



if we choose ul = i(l + ^) and vj = i(-l + 
The dispersion relation is given by 



Q{k) = JAl-ABl 



(A20) 



The above derivation and results apply exactly to the 
symmetric phase in the competitive case, only with the 



replacement ( | Jj_ | , | | ) 



'J±,—J'j_). The dispersion 



(|A20|I has the gapless point ci)(0,0, 0) in the AF phase 
where h — 0. It corresponds to the NG mode due to 
the spontaneous breaking of the U(l) symmetry. Hence 
Ah = ^(O, 0,0) can be regarded as the /i-induced gap. 
At the small field h, the angle variation around /i = 0: 
59 = sin~^(^^) — ^ci is estimated approximately as 



H 



(45'J)2 - iJ2 



Therefore expanding Af^{h,66) around h 
sees that the gap grows as 



(A21) 



0, one 



A, 



V45J 



\4SJ J 



1/4 



— y 

AS J 



1/2 



(A22) 

In contrast to the SSB phase (in the competitive case), 
the /i- induced gap opens as oc h^^^. Furthermore, 
it is remarkable that Eq. ljA22|l has no dependence on 
the inter-chain interactions. In fact, it is identical to 
the ID result [Eq. (3.17) in Ref.|l|. This is a reflection 
of the smoothness which was discussed in the final of 
Sec, mil At H = 0, we can obtain the simple exact result 
A,, = VISJh[l + h/{4:SJ)]^/^ from the dispersion HKM . 
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